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Abstract 

We study the radiative decay (f) — > 7r^7r'^7 within the framework of a phe- 
nomenological approach in which the contributions of cj-meson, p-meson and 
/o-meson are considered. We analyze the interference effects between different 
contributions and utilizing the experimental branching ratio for this decay we 
estimate the coupling constant gcji^'y 
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I. INTRODUCTION 



The radiative decays of mesons are valuable sources of information in low-energy 
hadron physics in areas such as quark model, SU(3) symmetry and the Okubo-Zweig-Iizuka 
(OZI) rule. In particular, radiative meson decays ir^TT^'y and — 71^777 can provide 
insight on the structure and properties of low-mass scalar resonances, since these decays 
primarily proceed through processes involving scalar resonances such as ^ /o (980)7 
ao (980)7, with the subsequent decays into and TT^rj'-f PIH- On the experimental 

side, the Novosibirsk SND |Q and CMD-2 |^ collaborations recently have reported very 
accurate results on these decays with the following branching ratios BR(0 — > 7r°7r°7) = 
(1.221 ± 0.098 ± 0.061) x 10'^ BR(0 ^ ttV) = (0-88 ± 0.14 ± 0.09) x 10"^ §, and 
BR(0 7r°7r°7) = (0.92 ± 0.08 ± 0.06) x 10~^ BR(0 ^ tt^) = (0-9 ± 0.24 ± 0.1) x 10"^ 
Q, where the first error is statistical and the second one is systematic. 

The low-mass scalar mesons with vacuum quantum numbers have fundamental impor- 
tance in understanding low energy QCD phenomenology and the symmetry breaking mech- 
anisms in QCD. In addition to /o(980) and ao(980), the existence of the scalar-isoscalar a 
meson as a broad tttt resonance, which has long been controversial, seems to be established. 
Increasing number of theoretical and experimental analyzes find a a-pole position near (500- 
1250) MeV A direct experimental evidence seems to emerge from an^ Sir 

decay channel observed by the Fermilab (E791) collaboration, where a meson is seen as a 
clear dominant peak with Mo-=478 MeV and ro-=324 MeV [0]. On the other hand, the 
nature and the quark substructure of these scalar mesons have not been established yet, 
whether they are conventional qq states tttt in case of a and KK in case of /o and 
ao [|10| molecules, or exotic multiquark g^g^ states [p!I] , p!2| have been a subject of debate. It 
has been noted that the radiative decay of meson to the scalar mesons can differentiate 
among various models of their structure IQ. 

The radiative decay process —>■ Tf^Tc^'j among other radiative decay processes of the type 

— ^ P^P^'-^ where V and P belong to the lowest multiplets of vector (V) and pseudoscalar 



1 



(P) mesons was studied by Fajfer and Oakes [0 using a low energy effective Lagrangian 



approach with gauged Wess-Zumino terms. They considered the contributions in which the 
virtual vector meson states dominate, and obtained the branching ratio for this decay as 
BR{(f) 7r°7r°7) = 3.46 x 10~^. The contribution of intermediate vector mesons to the 
decays — » P^P^'y was later considered by Bramon et al. []T1| using standard Lagrangians 
obeying SU(3)-symmetry, and they obtained the result BR{(f) 7r°7r°7) = 1.2 x 10^^ for the 
branching ratio of the decay. Bramon et al. later considered these decays 



within the framework of chiral effective Lagrangians, and using chiral perturbation theory 
they calculated the branching ratios for various decays of the type —>■ P^P^'-f at the 
one-loop level. They showed that the one-loop contributions are finite and to this order no 
counterterms are required. In this approach the decay 7r°7r°7 proceeds through charged- 
kaon loop and they obtained the contribution of charged-kaon loops to this decay as r(0 — >■ 
7r'^7r°7)=224 eV which is much larger than the intermediate vector meson state (VMD) 
contribution due to OZI rule. Considering the amplitudes of both VMD and kaon-loop 
contributions they obtained for the decay rate the value r(0 ^ 7r'^7r'^7)=269 eV, moreover 
they noted that OZI allowed kaon-loops dominate both the photonic spectrum and the decay 
rate. Radiative decays of meson were also investigated by Marco et al. [|r^ employing the 



techniques of chiral unitary theory developed earlier by Oiler |[T^. Using a chiral unitary 



approach, they included the final state interactions of two pions by summing the kaon-loops 
through Bethe-Salpeter equation. They obtained the branching ratio for the decay — >■ 
7r'^7r'^7 as BR{(f) n^'K^'-f) = 0.8 x 10^^. Moreover, they obtained the photon distribution 
as a function of the invariant mass of the vr^vr^ system and compared it with SND data. 
They noted that the shape of the experimental spectrum is relatively well reproduced with 
the — > /o7 contribution since /o meson is the important scalar resonance appearing in 
K^K^ Ti^-K^ amplitude. However, they also noted an appreciable strength for a possible 
(j) a'-f contribution in the spectrum. 

In this work, we follow a phenomenological approach and study the radiative — » TT^ir^'y 
decay by considering p-pole vector meson dominance amplitude as well as scalar cr-pole and 



/o-pole amplitudes. By employing the experimental value for this decay rate, we estimate 
the coupling constant g^a^-f- This coupling constant is an important physical input for studies 
of 0-meson photoproduction experiments on nucleons near threshold ||18|| . 



II. FORMALISM 

In our analysis of the radiative decay 7r*^7r°7, we proceed within a phenomenological 
framework and we do not make any assumption about the structure of /o meson. We 
note that the and /o mesons both couple strongly to K^K~ system, and therefore in 
our phenomenological approach we describe the (pKK- and /oi^TK-vertices by the effective 
Lagrangians 

CJkK = '9^KK<P^{K+d,K- ~K-d,K+) , (1) 

and 

Cf:^K = 9hKKMf,K+K~h , (2) 

respectively, which also serve to define the coupling constants g^ftRK and gj^KK- The effective 
Lagrangian for the (pKK-Yeiie'x is the one that results from the standard chiral Lagrangians 
in lowest order of chiral perturbation theory [|19|. The decay rate for the cj) K'^K'^ decay 



resulting from this Lagrangian is 

3/2 

(3) 
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Utilizing the experimental value for the branching ratio BR{(f) K^K~) = (0.492 ± 



0. 007) for the decay (f) K^K^ ||20|, we determine the coupling constant g^ftKK as g^KK = 
(4.59 ± 0.05). Furthermore, as a result of the strong coupling of both cj) and /o to K^K~ , 
independent of the nature and dynamical structure of /o, there is an amplitude for the decay 
~^ /o7 to proceed through the charged kaon-loop which we show diagrammatically in Fig. 

1, where the last diagram assures gauge invariance The amplitude of the radiative 



decay (p ~^ fol that follows from the diagrams in Fig. 1 is 



M{^^ /o7) = W^e^ ip.k, - g,.p . k) ' 9,KK^^f^KM,) ^^^^ 



(4) 



where (m, p) and (e, k) are the polarizations and four-momenta of the meson and the photon 
respectively, and a = M^/Mj^, b = Mj^jM"^. The function l(a,b) has been calculated in 
different contexts P, P^p2| , and it is given as 



I{a,b) 



2{a -h) (a - hf 



f 



f 



+ 



{a - by 



where 



arcsm 



ln(!^) - in 



X > 



(4x — 1)2 arcsin(2^ 



2v^^ ' ^ > 4 



i(l-4x)^[ln(^; 



X < 



(5) 
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The decay width for the radiative decay 
tude M.{(f) /07) as 



(6) 



/o7 can then be obtained from the ampli- 
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from which by using the experimental value BR{(j) /07) = (3.4 ± 0.4) x 10"'^ 
obtain the coupling constant QfoKK as g/t^KK = (4.13 ± 1.42). 

In our phenomenological approach, we assume that the radiative decay — s 
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(7) 



we 



proceeds through the reactions 



7r°7r°7, 



(77 — > vr vr 7, and — * /o7 



TT 7r 7. We write the total amplitude as the sum of the amplitudes of each reaction and 
this way we take the interference between different reactions into account. In order to 
describe the reaction — > /07 — > tt^tt^'j we again note that both and /q couple strongly to 
K^K~ , furthermore /o also couples strongly to Therefore, we assume that this reaction 

proceeds by a two-step mechanism with /q coupling to with intermediate KK states. We 
depict the processes contributing to the TT^ir^'-f decay amplitude diagrammatically in 



Fig. 2, where the last diagram in Fig. 2(c) results from the minimal coupling for gauge 
invariance. 

We describe the 0cT7-vertex by the effective Lagrangian 



nmt. _ 



cr 



(8) 



which also defines the coupling constant g,j,a'y. For the cttttt- vertex we use the effective 
Lagrangian 

1 



(9) 



The decay width of the a-meson that follows from this Lagrangian is given as 



= r(cr nn) 
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2M^ 



21 1/2 



(10) 



47r 8 

In our calculation, for given values of Ma- and Fo- we determine the coupling constant gcnrn 
by using the expression for T^. The (ppTc-vertex in Fig. 2(b) is conventionally described by 
the effective Lagrangian 



(11) 



The couphng constant g^p-,, is determined by N. N. Achasov and V. V. Gubin ||2^ as t^^pTr 
(0.811±0.081) GeV~^ using the data on the decay 
in Fig. 2(b) is described by the effective Lagrangian 



pvr TT+TT 71^ |2^. The p7r7- vertex 



(12) 



The coupling constant Qfy,^^ is then obtained from the experimental partial width of the 
radiative decay p 7r°7 as gpT,j = 0.69 ± 0.35. Finally, we describe the /o7r°7r°-vertex 
by the effective Lagrangian 

1 



(13) 



and the decay width Tf for the decay /o tttt that results from this effective Lagrangian 
is given by a similar expression as for Fo-. For a given value of Fj we use this expression 



to determine the coupling constant gfonn- Furthermore, in our calculation of invariant am- 
phtudes we make the replacement M — ^ M — |ir in /q-, a-, and p-propagators and use 
the experimental values Tp = (150.2 ± 0.8) MeV ^ for p-meson, and = (324 ± 21) 
MeV 1^ for cr-meson. However, the mass Mjg=980 MeV of /o-meson is very close to K^K~ 
threshold, and this induces a strong energy dependence on the width V of /o-meson. In 
order to take that into account different expressions for T can be used. A first possibility 
is to use an energy dependent width for /o 

r/o(g') = Tt{q') e{^^-2M^) + rj-(g2) e{^^-2MK) (i4) 

where is the four-momentum square of the virtual /o-meson. In this expression the width 
rjfo^(g^) is given as 



7o 



^ " (15) 
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and r'^-j-,(g^) with a similar expression. We use the experimental value r{'J^=40-100 MeV 



and we calculate T-^— from the effective Lagrangian given in Eq. (2). Another and 



widely accepted possibility is known from the work of Flatte [^. This amounts to extend- 



ing the expression for T^°-fr:{q^) below the KK threshold where J q^ — 4M|^ is replaced by 



KK 



isJAM]^ — q^ and thus V^^—[q^) becomes purely imaginary. In our work, we consider both 
possibihties. 

In terms of the invariant amplitude A^(E^, Ei)= M.a+M.h+M.c where M.a, M.b, and M.^ 
are the invariant amplitudes resulting from the diagrams a, b, and c in Fig. 2 respectively, 
the differential decay probability for -k^-k^'-) for an unpolarized 0-meson at rest is then 
given as 

dV 1 1 



M \\ (16) 



dE^dEi (27r)3 8M^ 

where E^ and Ei are the photon and pion energies respectively. We perform an average over 
the spin states of 0-meson and a sum over the polarization states of the photon. The decay 
width r(0 ^ 7r'^7r°7) is then obtained by integration 



1 C .max . C ^l.max. rJT^ 

T = =- dEj dEi-f—r (17) 

where now the factor (i) is included because of the 7r°7r° in the final state. The minimum 
photon energy is E^^min. = and the maximum photon energy is given as E^^max. = {M^ — 
4M|)/2M(^=474 MeV. The maximum and minimum values for pion energy Ei are given by 

1 



2{2E^M^ - Ml 



{-2E^M^ + 3E^Ml-Ml 



±EJ{-2E,M^ + Ml){-2E,M^ + - 4M2) } . 



III. RESULTS AND DISCUSSION 



The experimental full width of /o is F/g =40-100 MeV [^]. Since the coupling constant 
Qfj^^ry depends on the value of the width of /o, in order to estimate its effect on the coupling 
constant we take this width as = (70 ±30) MeV and Mjq=980 MeV in our calculation. If 
we use the form for T^^—{q^) that was proposed by Platte we are not able to reproduce 
the form of the invariant mass spectrum for the (j) tt^ti^'j decay. In this case the enhance- 
ment in the invariant mass spectrum appears in the central part rather than around the /o 
pole. A similar problem was also encountered by Bramon et al. in their study of the 
role of ao(980) exchange in the (p n^rj'-f decay. Therefore, in the analysis presented below, 
for rjp(g^) we use the form given in Eq. (14). Through the decay rate that results from the 
Lagrangian given in Eq. (13) describing /ovrvr-vertex we obtain the coupling constant gfoTrir 
as (7/Q7r7r=l-58. In order to determine the coupling constant g^a-y, we follow a similar proce- 
dure that we used in our previous works where we estimated the coupling constants gp^^ ||2^ 



and g^cr'y H^l- We use the experimental value of the branching ratio for the radiative decay 
— s> 7r'^7r'^7 in our calculation for this decay rate, and this way we arrive at a quadric 
equation for the coupling constant (^.^o-^, the coefficient of the quadratic term resulting from 
a meson contribution shown in Fig. 2(a), and the coefficient of the linear term from the 
interference of cr-meson amplitude with p-meson and /g-meson amplitudes shown in Figs. 
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2(b) and (c) respectively. Therefore our analysis produces for a set of values of a-meson 
parameters M„ and two values for the coupling constant g^a-^- We choose for the a- meson 
parameters the values M„ = (478 ± 17) MeV and T„ = (324 ± 21) MeV as suggested by 
the recent Fermilab (E791) experiment resulting for the coupling constant gcnm in the 
value (yfo-TTTT = 5.25 ± 0.32 using the expression for the decay rate T{a rni) given in Eq. 
(10). This way we obtain for the coupling constant g^pa-y the values g^pa-y = 0.064 ± 0.008 
and g^rj'y = 0.025 ± 0.009. We then study the invariant mass distribution for the reaction 
In Fig. 3 we plot the distribution dBR/dMT^T^ = {Mt^^^ / M^)dBR/ dE^ for the ra- 
diative decay 7r°7r°7 in our phenomenological approach choosing the coupling constant 
5'(/.o-7=0.064 and in Fig. 4 we show the same distribution for (yf0o-7=O.O25 as a function of the 
invariant mass M^^ of vr'^vr'^ system. In these figures we also indicate the contributions coming 
from the different reactions — * a7 — > 7r'^7r'^7, p^vr*^ — > 7r°7r°7, and /o7 — > 7r°7r°7 as 
well as the contribution of the total amplitude which includes the interference terms as well. 
The difference between the two total contribution curves is mainly in the interference region 
M^^ < 0.7 GeV above which /o-amplitude dominates the spectrum. From the analysis of 
the spectrum obtained with the coupling constants (7(^0-7=0.064 and g^„^=Q.Q2b in Figs. 3 
and 4, respectively, we may decide in favor of the latter value of the coupling constant g^pa-y 
and we may state that the experimental data within the framework of our phenomenological 
approach suggest that g^„^ = 0.025 ±0.009. In Fig. 5, we show the contributions of different 
amplitudes and the contributions of the interference terms in interference region M^^ < 0.7 
GeV for ^^,,^=0.025. 

On the other hand, the photoproduction of p°-meson on proton targets near threshold 
can be described at low momentum transfers by single one-meson exchange model. Friman 
and Soyeur ||28| showed that in this picture the p^'-meson photoproduction cross-section on 
protons is given mainly by a-exchange. They calculated pcr7- vertex assuming vector meson 
dominance of the electromagnetic current, and their result when described using an effective 
Lagrangian for the p(T7- vertex gives the value gp^^^y ^2.71 for this coupling constant. Later, 
Titov et al. in their study of the structure of 0-meson photoproduction amplitude based 
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on one-meson exchange and Pomeron-exchange mechanism used the couphng constant (^0^7 
which they calculated from the above value of gp^-y invoking unitary symmetry arguments 
as g^a'y ~0.047. Our result for this coupling constant (7^0-7 = 0.025 ± 0.009 is not in total 
agreement with this value. We note that in order to derive their result Titov et al. assumed 
that 0", /o, and are members of a unitary nonet. However, assignments of scalar states 
into various unitary nanets are not without problems and other possible assignments than 
used by Titov et al. have also been suggested In our work we do not make any 

assumption about the quark substructure of a and /o mesons and describe their couplings 
in a phenomenological framework. 

N. N. Achasov and V. V. Gubin performed a fit to the experimental data and they 
obtained the following values: the branching ratio with interference BR{(f) — > (/o7 + 7r°p°) — >■ 
7r°7r°7) = 1.26 X 10"^, the branching ratio of the signal BR{(f) — > /07 7r°7r°7) = 1.01 x 10""^, 
the branching ratio of the background BR{(f) p^n^ 7r°7r°7) = 0.18 x 10^^. If we use the 
coupling constant (7(^0-^=0.025 and M^=4:78 MeV ro-=324 MeV, we obtain for the branching 
ratios for different contributing reactions to the radiative decay 7r'^7r°7 the values 
/o7 7r°7r07) = 1.29 x 10-^ BR{^ ^ Tr^vr^) = 0.04 x 10'^, BR{(f) 

^ 7rV7) = 0.14xlO^^ BR{(j) (/o7+7rV) ^ 71^7) = 1.34xl0-^ BR{(f) (/07+ 
(T7) — » 7r°7r°7) = 1.16 x 10~^ and for the total interference term i?i?(interference)=— 0.25 x 
10~^. Our results are in reasonable agreement with those obtained in the analysis of N. N. 
Achasov and V. V. Gubin [|23|, with the difference that our results include contributions 
coming from the cr-pole amplitude as well as and /o-pole amplitudes. 

As we already noted, the spectrum for the decay 7r°7r°7 is dominated by the /q- 
amplitude, and the contribution coming from the cr-amplitude can only appreciably be 
noticed in the region M^-tt < 0.7 GeV through interference effects. Our analysis of these 
interference effects and our calculation for the decay rate of the radiative decay (f) 
suggests the value of the coupling constant g^j^a-y as g^j^a-y = 0.025 ± 0.009. 
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Figure Captions: 

Figure 1: Diagrams for the decay /07 
Figure 2: Diagrams for the decay (j) — > 7r°7r°7 

Figure 3: The 7r°7r° invariant mass spectrum for the decay (f) 7r°7r°7 for g<^o.^=0.064. The 
contributions of different terms are indicated. 

Figure 4: The n^n^ invariant mass spectrum for the decay (f) — > 7r°7r°7 for g<^o.^=0.025. The 
contributions of different terms are indicated. 

Figure 5: The contributions of p- and a- amphtudes and interference terms to the invariant 

mass spectrum of the decay (p — > n^n^j for g<^o-7=0.025. 
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